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There has been renewed interest in the frequency dependent specific heat of supercooled liquids 
in recent years with computer simulation studies exploring the whole frequency range of relaxation. 
The simulation studies can thus supplement the existing experimental results to provide an insight 
into the energy landscape dynamics. We here investigate a kinetic model of cooperative dynamics 
within the landscape paradigm for the dynamic heat capacity C(u>,T) behavior. In this picture, 
the /3-process is modeled as a thermally activated event in a two-level system and the a-process is 
described as a /3-relaxation mediated cooperative transition in a double well. The model resembles 
a landscape picture, apparently first conceived by Stillinger [Science 267, 1935 (1995)], where an 
a-process is assumed to involve a concerted series of /3-processes. The model provides a description 
of the activated hopping in the energy landscape in close relation with the cooperative nature of the 
hopping event. For suitable choice of parameters, the model predicts a frequency dependent heat 
capacity that reflects the two-step relaxation behavior. Although experimentally obtained specific 
heat spectra of supercooled liquids till date could not capture the two-step relaxation behavior, 
this has been observed in a computer simulation study by Scheidler et. al. [Phys. Rev. B 63, 
104204 (2001)]. The separation between the two peaks grows as the temperature drops, indicating 
the stringent constraint on the a-process due to the cooperativity requirement. The temperature 
dependence of the position of the low-frequency peak, due to the a-relaxation, shows a non-Arrhenius 
behavior as observed experimentally by Birge and Nagel [Phys. Rev. Lett. 54, 2674 (1985)]. The 
shape of the a-peak is, however, found to be temperature independent, in agreement with the 
simulation result. The high-frequency peak appears with considerably larger amplitude than the 
a-peak. We attempt a plausible reason for this observation that is in contrast with the general 
feature revealed by the dielectric spectroscopy. The relative amplitudes of the /3- and a-peaks in the 
present framework are found to depend on several characteristic features of the energy landscape, 
including the extent of cooperativity requirement for the a-relaxation and the asymmetry of the 
double well. 



I. INTRODUCTION 

Understanding the complex relaxation phenomena in 
supercooled liquids has motivated much scientific efforts 
over decades |l|, U, |3j . The measurement of frequency de- 
pendent specific heat, pioneered independently by Birge 
and Nagel and by Christensen in the year 1985 0, Hfl, 
opened up another approach to this goal. Zwanzig sub- 
sequently showed on the basis of linearized hydrodynam- 
ics that Cp(u>), the frequency dependent specific heat at 
constant pressure, could be directly related to the fre- 
quency dependent longitudinal viscosity ni{uj) 6J. The 
frequency dependent specific heat has, however, contin- 



ued to get explored from both experimental 0, H, @ and 
theoretical 0, O, 0, 0, Q] perspectives in the hope 
that specific heat spectroscopy would provide an insight 
into the energy landscape dynamics. 



In practice, one measures the frequency dependent spe- 
cific heat in the linear response regime following an arbi- 
trary small thermal perturbation that takes the system 
slightly away from the equilibrium c p (uj) is a linear 
susceptibility describing the response of the system to 
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this perturbation. One can, however, calculate the fre- 
quency dependent specific heat in terms of equilibrium 
fluctuation of energy following a relevant fluctuation- 
dissipation theorem derived explicitly by Nielsen and 
Dyre for a system whose dynamics is described by a mas- 
ter equation [l2j . In this work, we do so for a kinetic 
model of glassy dynamics that invokes the concept /?- 
organized-cc-process within the landscape paradigm. 

The measurements of frequency dependent specific 
heat currently suffer from a limitation that frequency 
range up to 10 4 Hz can only be probed with the presently 
available experimental setup. Therefore, the experimen- 
tally obtained frequency spectra capture only the a- 
relaxation regime of supercooled liquids. However, Scei- 
dler et. al. have recently carried out a computer simu- 
lation study of a system that models amorphous silica, 
where they could scan the whole frequency range of in- 
terest, revealing the two-peak structure with a notably 
dominant high-frequency peak [T^j . The high-frequency 
peak shows only a weak temperature dependence and 
has been ascribed to the vibrational excitations of the 
system. In their work ^4|, they have applied the Mori- 
Zwanzig projection operator formalism and made use of 
an exact transformation formula, due to Lebowitz et. al. 
|lfij |. to derive a relation between the frequency depen- 
dent specific heat c v (u) and the autocorrelation func- 
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tion of the temperature fluctuations in the microcanon- 
ical ensemble. This relationship, which is identical to 
the one derived independently by Nielsen in terms of a 
fluctuation-dissipation theorem, has allowed the determi- 
nation of c v (u>) from computer simulations in equilibrium 
[Tflj . A mode coupling theory (MCT) based calculation 
has also shown the two-step relaxation behavior in the 
predicted frequency spectrum of the specific heat |17| . 
The dominance of the high-frequency peak is, however, 
not evident in this work. 

The measurements of frequency dependent specific 
heat in specific heat spectroscopy supplemented by com- 
puter simulation studies may prove to be useful in pro- 
viding insight into the landscape dynamics of supercooled 
liquids. The landscape paradigm has been widely used to 
elucidate d yna mics of liquids in the supercooled regime 
\}M fl9t Em l2lj . This framework involves the division of 
the multidimensional configuration space into so called 
metabasins on the basis of a transition free-energy cri- 
terion. Two vastly different timescales thus get entailed, 
the smaller one due to motions within the metabasins and 
the longer one due to exchange between the metabasins 
involving much larger free-energy of activation. In par- 
ticular, the /3-processes are visualized to originate from 
activated dynamics within a metabasin, while escape from 
one metabasin to another is taken to describe an a- 
process |2jj. See Figure 1 for a schematic represen- 
tation of the two processes. It is important to note here 
that the breakdown of the MCT ElliJ] is ascribed 
to the dominance of relaxation by these thermally acti- 
vated hopping events |25j, which are unaccounted for in 
the ideal version of MCT. Recent computer simulation 
studies have further revealed that hopping is a highly 
cooperative phenomenon promoted by many body fluc- 
tuations [2f| H3, Hlj; hopping of a tagged particle 
is often preceded by somewhat larger than normal, but 
still small amplitude motion of several of its neighbors 
[29| . A rather different stringlike cooperative motion has 
also been found to occur in a model glass-forming liquid 
|3C|. In a very recent simulation study that investigates 
the role of cooperativity in reorientational and structural 
relaxation in a supercooled liquid both on the potential 
energy landscape and in the real space, it has been fur- 
ther proposed that relaxation involves activation of the 
system to the complex multidimensional region on the 
borders of the basins of attraction of the minima for an 
extended period |3lj . 

In the present work, we employ a kinetic model of 
glassy dynamics that attempts to provide a description 
of the activated hopping within the landscape paradigm 
in close connection with the cooperative nature of the 
hopping event. We follow a procedure, as outlined by 
Nielsen and Dyre [l2|. to compute the frequency depen- 
dent heat capacity C(v,T) for our model system. For 
suitable choice of parameters, the model predicts a fre- 
quency dependent heat capacity that captures many of 
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FIG. 1: A schematic representation of the potential energy 
landscape showing motions within and between metabasins. 



the features of the two-step relaxation behavior in su- 
percooled liquids. The predictions of our model are in 
good qualitative agreement with the available experimen- 
tal and computer simulation results. 

The outline of the paper is as follows. In the next 
section we describe the model. Section III provides the 
theoretical treatment. We present the results along with 
discussion in section IV. Section V concludes with a sum- 
mary of the results and a few comments. 

II. DESCRIPTION OF MODEL 

We model a /3-process as an activated event in a two- 
level system (TLS). We label the ground level of a TLS 
as and the excited level as 1. The waiting time before 
a transition can occur from the level i{= 0, 1) is assumed 
to be random and is given by the Poissonian probability 
density function: 



ipi(t) = —exp(-t/n), 

Ti 



i = 0,l, 



(1) 



where Ti is the average time of stay at the level i. If 
Pi(T) denotes the canonical equilibrium probability of 
the level i of a TLS being occupied at temperature T, 
the equilibrium constant K(T) for the population in two 
levels at temperature T is given by the following relation 
that obeys the detailed balance: 



K(T) 



pi(T) = ri(T) 
Po(T) r (T) 



exp[-e/(k B T)], (2) 



where e is the energy separation between the two levels 
in a TLS, and ks is the Boltzmann constant. Here the 
level is taken to have a zero energy. 

Within the framework of the present model, a 
metabasin is characterized by an Np number of such non- 
interacting two-level systems (TLSs). A given minimum 
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/ displacements are required simultaneously among its 
neighbors. The activated dynamics within a TLS may 
well represent small rotations ^j- ^ n the case of poly- 
mer melts which exhibit glassy behavior, the /^-relaxation 
may involve the motion of side chains. This picture ap- 
parently differs from the one drawn by Dyre |32j| , who has 
argued that large-angle rotations are " causes" and small- 
angle rotations are "effects". The present picture, how- 
ever, contains Dyre's one in the sense that small-angle 
rotations indeed occur following a large-scale jump mo- 
tion for the completion of relaxation as evident in Figure 
2. The present model is built on a rather symmetrical 
picture that also necessitates small-angle rotations for a 
large-angle rotation to occur. 



FIG. 2: A schematic representation of the model under con- 
sideration. The horizontal lines within a well represent differ- 
ent excitation levels. Note that the energy levels are in general 
degenerate, as they correspond to the sum of the energies of 
individual TLSs in the collection. 



number among the total number Np of TLSs must si- 
multaneously be in the excited levels for the occurrence 
of an a-process. We here concentrate on two adjacent 
metabasins, which we label as 1 and 2 and together 
call a double well. Figure 2 shows a schematic dia- 
gram of two adjacent metabasins with illustration of dy- 
namics within and between them. The respective num- 
bers of TLSs that comprise the metabasins are Ni 1] and 



N 



(2) 







For a collection of N { 1 ' (i = 1,2) TLSs, a vari- 
able Q(t), (j = 1,2, ,N^) is defined, which takes 

on a value if at the given instant of time t the level 
of the TLS j is occupied and 1 if otherwise. Q(t) 
is thus an occupation variable. The collective variables 
Qi(t) (z = 1, 2) are then defined as 



i(t) 



Em 



(3) 



Qi (t) is therefore a stochastic variable in the discrete in- 
teger space [0,Ng']. Qi(t) serves as an order parameter 
for dynamical change involving metabasin i. Here an a- 
process is assumed to occur only when all the (3-processes 
(TLSs) in a metabasin are simultaneously excited, i.e., 



when Qi 



N 



(i) 




There is a finite rate of transition k 



from each of the metabasins when this condition is sat- 
isfied. Within the general framework of the model, the 
double well becomes asymmetric when ^ A/1 2 , as 
shown in Figure 2. 

It is worth while to note the correspondence of the 
present description with real physical processes occurring 
in glass-formers. The a-process may correspond to large- 
scale hopping of a particle. For this hopping to occur, 
however, many small reorientations / rearrangements 



III. THEORETICAL TREATMENT 

From a theoretical point of view, the treatment of fre- 
quency dependent heat capacity can be carried out by 
employing the linear response assumption. Following 
Nielsen and Dyre 0, the frequency dependent heat ca- 
pacity C(uj,T) of our system at temperature T can be 
given by 



CV,T) 



<E 2 (T) > 



k B T 2 



k B T 2 



dt 



<E(t,T)E(0,T)>,(4) 



where s = ioj, ui being the frequency of the small oscil- 
lating perturbation, i = \/—\, and the angular brack- 
ets denote an equilibrium ensemble averaging. E(t,T) 
stands for the total energy of the system at time t and 
temperature T and is given by 



A' 



E(t, t) = *. T ) i N f ] - N ? + ») e 



N\. 



^P 2 (n-t,T) n e. 



(5) 



n=0 



Here the lowest level of the well 2 is taken to have zero 
energy and Pi(n;t,T) denotes the probability that the 
stochastic variable Qi takes on a value n in the i-th well 
at time t and temperature T. The evolution of these 
probabilities obeys the master equation |33l |: 



dPj(n;t,T) 
dt 



= [(JV«-n + l)/7 ti (T)]P i (n-l;t,T) 

+ [(n + l)/ri(T)]Pi(n+l;t,T) 

- [(JVf-n)/r (r)]Pi(n;t,T) 

- (n/Ti(T))P i (n;t,T)-kS j 



,1V 



+ k 6 (i ±i) 5j^±i Pj(n;t,T), 



P(n;t,T) 
(6) 



where the '+' and '— ' signs in the indices of the Kronecker 
delta are for i = 1 and 2, respectively. 
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One can have the following compact representation of 
the set of equations given by Eq. © for all possible n and 
i values 



dP{t,T) 
di 



A(T)P(i,T), 



(7) 



where Pi(n; t, T) for n = 0, 1, A^ 13 and P 2 (n; t, T) for 

n = 0, 1, ...,Ng ' together comprise the elements of the 
column vector P(£, T) and A is the transition matrix of 

order N = Nf ] + + 2. If G T (i, t\j, 0) be the Green's 
function that gives the probability to be in the state i 
at a later time t given that the system is in the state 
j at time t' = 0, the temperature being kept constant 
at temperature T, the matrix of Green's functions also 
satisfies the rate equation 



dG T (t) 
~dt 



= A(T)G T (t) 



(8) 



with the initial condition Gy(Q) = I, where I is the iden- 
tity matrix of order N. In terms of Green's functions, 
one can then rewrite the energy autocorrelation function 
as 

N N 

< E{t,T)E{0,T) >=J2 J EGT(i,t\j,0)E i E j P eq (j,T), 

2=1 j = l 

where P eq (j, T) is the equilibrium probability of the state 
j at T. We write GT(i,s\j) as the Laplace transform of 
G T (i,t\j,0): 

OO 

G T (i,s\j) = J dte- st G T (i,t\j,0). (10) 
o 

The frequency dependent heat capacity is then given by 
< E 2 (T) > 



C(u,T) = 



N N 



k B T 2 



i=i j=i 



(ii) 



The computational procedure involves the numerical 
evaluation of the Green's functions by an inversion of 
matrix: 



G T (*) = (*I-A(T))- 



(12) 



and that of P eq (T) from the eigenvector corresponding 
to the zero eigenvalue of A(T). In the next section, we 
present the results with discussion and note the relevance 
of our results. 
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FIG. 3: Frequency dependence of the imaginary part of the 
dynamic heat capacity C"(uj*, T*) for our model system with 
A 7 ! = 3 and A^ 2 ' = 5, at three dimensionless reduced tem- 
peratures T* — 0.9 (solid line), T* = 0.8 (dashed line), and 
T* = 0.7 (dot-dashed line). Temperature T is scaled by 
the melting temperature T m to have a reduced temperature 
T* — T/Tm- Frequency is also scaled by the inverse of ri(TWi) 
to get a dimensionless reduced frequency lu* — u>Ti(T m ). As 
in the earlier work. [l5|l we set k = 2.0 in the reduced units, 
and e = 2ksT m and e\ = 18fesT m , the latter being the en- 
ergy barrier to the transition from the level 1 in a TLS. The 
presence of a fixed energy barrier to transition from one level 
to the another within a TLS is expected to impart an Ar- 
rhenius temperature dependence of the /3-relaxation within 
the present framework. The laboratory glass transition tem- 
perature T g occurs at a temperature around (2/3)T m .Q All 
temperatures investigated here therefore fall between T m and 
(2/3)T m . The same set of parameter values has been used 
for all calculations in the present work unless it is specifically 
mentioned otherwise. 



metabasin, the energy separation e between the two levels 
of a TLS, the energy of activation e* for barrier crossing 
within a TLS, the energy asymmetry A between the two 
adjacent metabasins, and the critical number N c of TLSs 
required to be in the excited levels at a particular time 
for the a-relaxation to occur. The choice of these pa- 
rameters has been kept simple in this work though at 
the expense of being ad hoc at least in some cases. For 
example, we have taken N c to be equal to Np and the 
value of e has been taken of the order of fcsT m , T m being 
the melting temperature. An approximate estimate of 
Np has been taken from simulations results |28l |29| . The 
value of the activation energy is a rather difficult guess. 
We have used the guidance provided by an earlier work 
of ours [15| . 

IV. RESULTS AND DISCUSSION 



The dynamical response of the system in the present 
framework is expected to be determined by a set of pa- 
rameters that includes the number Np of TLSs in a 



In Figure 3, we show the frequency dependence of 
the imaginary part of the heat capacity C" (oj* ,T*) cal- 
culated for our model system at three different temper- 



5 




-4 



3 -5 

DC 
O 



-6 




-0.5 1 — ■ 1 ■ 1 ■ — 1 — ■ — 1 — ■ — 1 — ■ — 1 — ■ — 1 

-10 -8 -6 -4 -2 2 4 

logoo 

FIG. 4: Frequency dependence of the real part of the dynamic 
heat capacity G'(u)*,T*) for our model system at three re- 
duced temperatures T* = 0.9 (solid line), T* = 0.8 (dashed 
line), and T* = 0.7 (dot-dashed line). 



atures. The two-peak structure corresponding to the bi- 
modal relaxation behavior as conceived in the model is 
evident at all temperatures investigated. The peak at 
high frequencies corresponds to the (^-relaxation and is 
remarkably dominant. The low-frequency peak is due to 
the a-relaxation that occurs on a longer timescale. Note 
that the separation between the position of the (3-peak and 
that of the cx-peak grows as temperature is lowered. This 
is indicative of the stringent constraint on the a-process 
that arises from the cooperativity requirement. 

The features of the two-step relaxation also get re- 
flected in the frequency spectrum of the real part of the 
heat capacity C'(ui*,T*) as shown in Figure 4 at three 
different temperatures. This is expected on the basis 
of the Kramers-Kronig relation that relates the real and 
imaginary parts. Whenever lo~ 1 is on the order of the 
timescale of a characteristic relaxation process, the sys- 
tem takes up energy inducing an increase in the real part 
of the heat capacity around that frequency. The domi- 
nance of /3-relaxation is again evident from a much larger 
increase at the high frequencies. The low-frequency limit 
corresponds to the static heat capacity of the system. 

Let us now discuss the temperature dependence of the 
positions of the peaks as they appear in the frequency 
spectrum of C"(u>*,T*). The temperature dependence 
of the a-peak position oj* a (T*) in the reduced scale 
is shown in Figure 5 in a log oj* a (T*) versus 1/T* 
plot. The data fit well to the Vogel-Fulcher-Tammann 
(VFT) equation: w* a (r*) = oJoy FT exp[-A VFT /{T* - 
To.vft))- The small curvature of the fitted curve is, 
however, notable. This has been observed experimen- 
tally as well Q. As in Ref. [4], we show in Figure 5 a 
second fit to the data with a scaling law: lj* a (T*) — 
^o,sci[{T* - T 0tSC i)/T 0iSC i]. The two fits are nearly in- 
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FIG. 5: The a-peak frequency w*„, on a logarithmic scale, 
versus the inverse temperature 1/T*. The solid line and 
dashed line correspond to two nearly indistinguishable fits 
to the data with the Vogel-Fulcher-Tammann equation and 
a three-parameter scaling law (see text). From the fits, 
To, vft = 0.154 and Tq.sc! = 0.475 in the reduced units. The 
inset shows the temperature dependence of the /3-peak fre- 
quency a in a logut a versus the 1/T* plot. The solid line 
is a linear fit to the data with a slope of 7.912 in the reduced 
temperature units. 



distinguishable. The inset of Figure 5 shows the tem- 
perature dependence of the /^-relaxation peak position 
w*p(T*) in a similar plot. The linear fit corresponds to 
the Arrhenius behavior as expected from the supposition 
of the model. 

Since the experiments till date could probe only the 
a-relaxation regime of the frequency spectrum, we have 
investigated the features of the a-relaxation peak as pre- 
dicted by our model system in a bit more detail. In Fig- 
ure 6, we concentrate on the a-relaxation regime of the 
frequency spectrum of C" (oj* , T*) at different tempera- 
tures to get an enlarged view of the a-relaxation peak. It 
is evident that the area under the a-peak becomes smaller 
as temperature drops. This implies that the cooperative 
component of the configurational part of the heat capac- 
ity diminishes with decreasing temperature. See Ref. [14] 
for a detailed discussion. Although the height and the po- 
sition of the a-peak have been found to change with tem- 
perature, its shape seems to show no such dependence. 
This can be evident if one constructs a master plot by 
scaling both the height of the peak and the frequency. 
This is demonstrated in the inset of Figure 6 where 
we plot C"(uj*,T*)/C"(oj* i01 ,T*) versus log(oj* /oj* a ) at 
different temperatures for the frequency range of inter- 
est. The curves at different temperatures appear to be 
superimposed on each other with negligible error. It is 
therefore reasonable to conclude that the shape of the a- 
peak here does remain independent of temperature. This 
is also in remarkable agreement with the observation in 
the computer simulation study reported in Ref. [14]. 



FIG. 6: The frequency spectra of the imaginary part of the 
dynamic heat capacity C" (lu* , T*) at different temperatures 
showing the low-frequency a-peak only. From the right to the 
left T* = 0.95,0.9,0.85,0.8,0.75,0.7, respectively. The inset 
shows the same curves, in dotted lines, scaled by the height 
of the respective peaks versus log(u* /ut a ). The solid line is 
a fit to the curve for T* =0.7 with a three-parameter equa- 
tion that is a good frequency domain representation of the 
time domain Kohlrausch- Williams- Watts stretched exponen- 
tial form. The dashed line is a fit to the same curve with the 
Debye response function. The curves along with the two fits 
are nearly indistinguishable. The details of the fit are given 
in the text. The range of frequency shown here is restricted 
at the high frequency side due to the presence of the /3-peak 
that leads to a break down of the fitting. 

The results presented here till now have been for a sin- 
gle set of Np~\ , e and e\ values. It is important 
to know how the predicted results are dependent on the 
choice of model parameters, in particular on ivi and 

(2) 

Nf', before one judges the merits of the model. To this 
end, we have explored the parameter dependence of the 
predictions of our model. The two-step relaxation behav- 
ior has been found to be revealed in the predicted heat 
capacity spectra of the model for a few other choices of 
{Nf\Nf ] } as well, specifically, for {4,6}, {5,7}, and 
{6,8}. For each of these parameter sets, a similar non- 
Arrhcnius behavior of the temperature dependence of the 
a-peak position has been observed as shown in Figure 7. 
The Vogel temperature Toy ft obtained from the VFT 
fit shifts to higher values along the series as one would 
expect from an increasingly stringent cooperativity re- 
quirement. Although the magnitude of the shift is small, 
the trend is quite clear. The shape of the a-peak has also 
been found to remain independent of temperature along 
the series (data not shown). For {6,9}, the amplitude 
of the a-peak gets diminished to a large extent by the 
time T* — 0.7 is reached, while the two-peak structure 
does not appear at all even at T* = 0.95 for {6, 10}. 
The two-step relaxation behavior also does not get evi- 
dent from the predicted dynamic heat capacity spectra 
in the temperature range explored here for even {2, 6}. It 



FIG. 7: The a-peak frequency w*„, on a logarithmic scale, 
versus the inverse temperature 1/T* for three different sets 
of {N^ , } values. For each set, the solid line and 
dashed line correspond to two nearly indistinguishable fits 
to the data with the Vogel- Fulcher-Tammann equation and a 
three-parameter scaling law. From the fits, Tq,vft = 0.174 
and T 0iSci = 0.485 for {4,6} (top), T 0tVFT = 0.197 and 
To.sd = 0.496 for {5,7} (middle), and T ,vft = 0.203 and 
T , sc i = 0.498 for {6, 8} (bottom). 



follows therefore that a difference of 4 in Ng and Ni 

with the set of e and e\ values used here is large enough 
to suppress any manifestation of the two-step relaxation 
behavior within the present framework. On considering 
varying asymmetry of the two adjacent metabasins, we 
further make an interesting observation as illustrated in 
Figure 8. Here, we have varied with Np , e and e\ 
held fixed. For a symmetric double well, only one peak 
appears in the frequency spectrum of the imaginary part 
of the dynamic heat capacity. The a-peak disappears as 
an a-process is inconsequential from an energy consider- 
ation for the symmetric double well. The amplitude of 
the a-peak is the largest with the least asymmetry and 
gradually diminishes with growing asymmetry before it 
gets suppressed completely. The results suggest that the 
relative amplitudes of the a-peak may provide an insight 
into the energy asymmetry between the metabasins in 
the energy landscape of the system. 

In the inset of Figure 6, we have also shown a fit to the 
curve for T* = 0.7 with a three-parameter equation that 
is a good frequency domain representation of the time do- 
main Kohlrausch- Williams- Watts (KWW) stretched ex- 
ponential form [34||. From the fit parameters we obtain 
the stretching parameter (3kww = 0.94. Such a /3kww 
value implies very weakly nonexponcntial behavior. In 
fact, a fit with the Debye response function that corre- 
sponds to a single exponential behavior is found to be rea- 
sonable and for some of the other sets of {ivi 1 , Np}, is 
found to be as good as the one with the three-parameter 
equation corresponding to a stretching parameter very 



log CO 

FIG. 8: The frequency spectra of the imaginary part of the 
dynamic heat capacity C" (lo* ,T*) at temperature T* = 0.95 
for four choices of with kept fixed: (a) {3, 3} (solid 
line), (b) {3,4} (dotted line), (c) {3,5} (dashed line), and (d) 
{3, 6} (dot-dashed line) . 
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FIG. 9: The frequency spectra of the imaginary part of the 
dynamic heat capacity C" (u>* ,T*) at temperature T* = 0.9 
for four sets of {N^,Nf } } values with Nj^ - held 
fixed: (a) {4,6} (dotted line) , (b) {3,5} (solid line), (c) {2,4} 
(dotted line), and (d) {1,3} (dot-dashed line). 



close to unity. This is reasonable as spatially hetero- 
geneous domains, which is believed to be the primary 
reason for the stretched exponential relaxation in super- 
cooled liquids j3^, |3(| |33, has not been considered in 
the present calculations. The heterogeneous dynamics in 
different domains can be included either through a dis- 
tribution of e (the separation between the energy levels 
within a TLS) or through a distribution of barrier height 
for transition from one level to the other within a TLS. 
When the heterogeneity is included, the exponent Pkww 
is expected to decrease considerably as indeed found in 
the treatment of structural relaxation within a similar 
model where we have considered a distribution of e |3Sj . 

The remarkably dominant /3-peak as predicted by our 
model merits further consideration. In order to trace 
back its origin, we note that the constraint of cooper- 
ativity on the a-relaxation allows the system to take 
up energy only through the localized /^-processes unless 
the condition is satisfied. The system takes up more 
energy through the excitation of individual /3-processes 
(that collectively bring about the a-process) than the a- 
process itself which corresponds to a transition from one 
metabasin to another. This argument is substantiated 
in Figure 9, where we focus on the amplitude of the 
/3-peak relative to that of the a-peak in the frequency 
spectrum of the imaginary part of the dynamic heat ca- 
pacity. Note that the a-peak is not distinctly observed 
for {1, 3} as the separation between the two peaks is not 
enough. A decrease in the requirement of cooperativity 
for the a-process is found to result in a fall in the relative 
amplitude of the /3-pcak. One can, therefore, define a pa- 
rameter Q?P as the ratio of the amplitude of the /3-peak 
to that of the a-peak, which is in general temperature 
dependent. That is, Q£f(T) = C"(w Ptf) ,T)/C"(u p , a ,T). 



The parameter Q£f(T) may serve as a measure of the 
cooperativity needed for a-relaxation to take place. 

Let us finally consider a limiting case where the /3- 
processes are not associated with any energy cost, i.e. 
e = 0. We then describe a /3-process as an activated 
event in a two-state system (rather then a two-level sys- 
tem) with a nonzero energy of activation e* for barrier 
crossing. We further assume a nonzero value of A and 
one of the fcj's to be temperature independent, where fcj 
is the rate of transition from the metabasin i (See Fig- 
ure 2). In particular, we assume &2 to be independent 
of temperature while k\ = ki exp[— A/(fesT)]. From the 
results we have already presented here, one would expect 
then the imaginary part of the dynamic heat capacity to 
show only one peak. This has indeed been observed as 
shown in Figure 10. 

V. CONCLUSION 

Let us first summarize the present work. We have em- 
ployed a kinetic model of glassy dynamics that considers 
cooperativity through the constraint of [3-organized-a- 
process within the landscape paradigm. The two-step 
relaxation behavior as conceived in the model gets re- 
vealed in the frequency dependent heat capacity of the 
model for a reasonable range of parameter values. The 
analysis of the predicted dynamic heat capacity spectra 
suggests the following emergent features of our model: 
(j) The a-peak frequency has a non-Arrhcnius temper- 
ature dependence, (ii) The shape of the a-peak as it 
appears in the frequency spectrum of the imaginary part 
of the dynamic heat capacity is invariant in temperature, 
(m) The amplitude of the /3-peak is considerably larger 
than that of the a-peak. (iv) In a limiting case where 
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FIG. 10: The frequency spectra of the imaginary part of the 
heat capacity C"(uj*,T*) at different temperatures for the 
limiting case considered with Nf' = 3, = 5, e = 0, e* = 
l8kBT m , and A = 4fcsT m . We set = 2.0 in the reduced 
units and the temperature dependence is kept contained in ki . 
From the right to the left T* = 0.95,0.9,0.85,0.8,0.75,0.7, 
respectively. 



the /3-process involves an activated transitions between 
two states of the same energy while the two adjacent 
metabasins have an asymmetry in energy, the frequency 
spectrum of the dynamic heat capacity exhibits only a 
one-peak (one-step) structure. 

While the first (the non-Arrhenius temperature depen- 
dence) of this lot conforms to a number of experimental 
results, the second one has been observed in a recent com- 
puter simulation study |l4j |. The third one, however, does 
invite a few comments. Although a remarkably domi- 
nant high-frequency peak has been reported in Ref.[14], 
it is rather ascribed to the vibrational excitations of the 
system. We, however, note that the predominance of the 
[3 -peak is not in contention of being contradicted with any 
published experimental data on the dynamic heat capac- 
ity. In this context, it may be reasserted that the /3- 
relaxation regime is still beyond the scope of the specific 
heat spectroscopy since the existing experimental setup 
can scan a frequency range up to only 10 Hz. Neverthe- 
less, such a dominant /3-peak is unlikely to be observed in 
typical fragile liquids where even for a deeply supercooled 
state, the signature of the /^-relaxation in the frequency 
spectrum of the dynamic structure factor is only a weak 
secondary peak at higher frequencies. 

The present work suggests that the relative amplitudes 
of the two peaks may provide insight into the microscopic 
mechanism of the relaxation processes in deeply super- 
cooled liquids. In view of this, further measurements of 
frequency dependent specific heat of various glass form- 
ing liquids will certainly be worthwhile. In a recent sim- 
ulation study of computer model of amorphous Ni^iBg, 



van Ee et. al. have shown that the hopping mode is not 
only collective but appears to involve rather large number 
of neighbors • This study seems to support the pic- 
ture that vibrational coherence among neighbors is a pre- 
requisite for large scale hopping. In such a scenario, our 
/3-process may be identified with such vibrational motion 
while hopping is the jump between the two minima shown 
in Figure 2. We further argue that a predominant /3- 
peak would be a characteristic feature of Stillinger's pic- 
ture which assumes an a-process to involve a concerted 
series of /3-processes |40|. However, proper characteriza- 
tion of the /3-processes remains unclear. It may be noted 
here in this respect that the activated dynamics with an 
Arrhenius type temperature dependence characterizing 
the /?-process in the present work is the typical feature 
of what is called the 'slow' /3-process |40j. However, the 
manifestation of this dynamical process with a predomi- 
nant peak in the frequency spectrum of the dynamic heat 
capacity is more likely to be relevant with an even faster 
/3-process typically appearing in the mode-coupling the- 
ory predictions. 

Several other comments on the present work are in or- 
der. First, the present model can be taken to belong to 
the class of kinetically constrained models that attempts 
to provide a description of glassy dynamics by imposing 
dynamical constraints on the allowed transitions between 
different configurations of the system, while maintaining 
the detailed balance |4lj . In particular, our model resem- 
bles the models of hierarchically constrained dynamics of 
glassy relaxation, due originally to Palmer et. al. [42| . 
in the spirit that brings in cooperativity. Second, the 
high-frequency peak for real liquids is likely to draw con- 
tribution from /3-relaxation as well as vibrational excita- 
tions. In some cases, specifically at low temperatures 
close to the glass transition, it is possible that these 
two can be sufficiently separated to give rise to an ad- 
ditional peak in the imaginary part of the frequency de- 
pendent specific heat. Third, it is imperative to com- 
pare the two-peak structure of the frequency spectrum 
of the imaginary part of the dynamic heat capacity with 
its dielectric analog. In contrary to the prediction our 
model makes on dynamic heat capacity behavior, it is 
rather the predominance of the a-peak that has been 
observed in a vast body of experimental data on dielec- 
tric relaxation. While the fluctuation in energy within 
the present framework of the model that has an energy 
landscape picture at the backdrop translates easily into 
the calculation of the frequency dependent heat capacity, 
the model as such does not allow us to calculate the fre- 
quency dependence of the dielectric constant. The lat- 
ter needs further development of the model. However, 
the well-known bimodal frequency dependence of the di- 
electric relaxation in supercooled liquids can be at least 
qualitatively understood from the present description of 
/3- and a-processes. We essentially follow the description 
of Lauritzen and Zwanzig in assuming that a /3-process 
can be taken to correspond to a two-site angular jump 
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of individual molecules by a small angle around some 
axis 01 . These individual, uncorrelated angular jumps 
lead to a partial relaxation of the total electric moment 
M(t) of the whole system (note that M(f) is the sum of 
the dipole moment of the individual molecules). The di- 
electric susceptibility spectrum can be obtained from the 
auto-time correlation function of M(t) by using the linear 
response theory |44|. Since M(i) is a sum of a relatively 
large number of individual dipole moments, the former 
is a Gaussian Markov process and thus the time corre- 
lation function of the /^-relaxation mediated part must 
decay exponentially. As noted earlier, this /3-relaxation 
mediated decay is incomplete because all the jumps are 
small and restricted. Thus, it is fair to assume the follow- 
ing form for the auto-time correlation function of M(t) 

C M (t) = < Mp > exp{-t/T P ) 

+ (< Ml >-< Mj >) exp(-t/T a ), (13) 

where rp and r Q are the time scales of /3-relaxation and a- 
relaxation, respectively. In the above equation < Ml > 
is the value by which the mean-square total dipole mo- 
ment decays due to (3 relaxation alone from the initial 
value of < M 2 >. The rest of the decay (i. e. from 
< M D 2 > — < M% >) to zero occurs via the a-process. 
This suggests that with well separated rp and r , one 
would observe bimodal dispersion. However, the calcu- 
lation of < Mp > would require a more detailed model 
than the one attempted here. 



The present study suggests several future problems. 
First, it would be interesting to investigate the frequency 
dependence of the specific heat of several molecular liq- 
uids of varying fragility in computer simulations. In some 
cases, one should be able to discern a three-peak struc- 
ture. Second, it would be also of interest to investigate 
if the frequency dependent specific heat can be used to 
study the so called Boson peak which has drawn much 
attention in recent times |45j. Third, a comparison of 
the frequency dependence of specific heat with that of 
the dynamic structure factor of certain model systems 
would provide us insight to decide on whether the micro- 
scopic mechanisms of relaxation vary for different modes. 
Finally, a generalization of the present model or an al- 
together different model to describe the frequency de- 
pendence of dynamic heat capacity and that of dielectric 
response in the same framework would be a worthwhile 
undertaking. 
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